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Abstract 

The masses of 1 Sq and 3 S\ mesons containing a single heavy quark are com- 
puted in the quenched approximation. The light quark action and gauge field 
action are both classically-improved and tadpole-improved, and the couplings 
to the heavy quark are organized by the 1/M expansion of tadpole-improved 
NRQCD. At each of two lattice spacings, near 0.22fm and 0.26fm, meson 
masses are obtained for heavy quarks spanning the region between charmed 
and bottom mesons. Results up to 0(1/M), 0(1/M 2 ) and 0(1/M 3 ) are dis- 
played separately, so that the convergence of the heavy quark expansion can 
be discussed. Also, the effect of each term in the 0(1/M 3 ) contribution is 
computed individually. For bottom mesons the 1/M-expansion appears to be 
satisfactory, but the situation for charmed mesons is less clear. 
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I. INTRODUCTION 



Long-distance, nonperturbative QCD interactions can be studied numerically by dis- 
cretizing space-time, if the lattice spacing "a" is sufficiently small to allow a matching to 
perturbative QCD. To include the effects of a quark whose inverse mass is smaller than the 
lattice spacing, it is natural to use an effective Lagrangian which is ordered by powers of 
the inverse quark mass. Two different notations which have been used for this expansion 
are heavy quark effective theory (HQET) || and nonrelativistic QCD (NRQCD) g,|§. 

In the present work, the NRQCD formalism is used to study the masses of the ground 
state ( 1 5'o) and first excited state ( 3 5i) heavy-light mesons, i.e. mesons containing a single 
heavy quark. These masses are well-known experimentally 0, and have been previously 
determined from lattice NRQCD [f^| up to 0(1/M 2 ). The primary goal of the present 
research is to extend the calculation to 0(1/M 3 ), and to display the effects of each new 
term individually. This provides an indication of the convergence of the 1/M expansion 
of lattice NRQCD. Of particular interest are the physically-relevant cases of charmed and 
bottom mesons, both of which will be discussed herein. 

To allow the use of coarsely-spaced lattices, the light-quark and gauge terms in the action 
will be classically-improved, and the entire action will be tadpole- improved. Not only does a 
larger lattice spacing imply speedier simulations, but it also means that the minimum heavy 
quark mass of NRQCD (which is of the order of the inverse lattice spacing) is reduced. This 
opens the possibility of using NRQCD to study charmed mesons. 



Lattice NRQCD has been used extensively for studies of quarkonium | ID| , where the 
1/M expansion is replaced by a velocity expansion, and it was concluded that the velocity 
expansion for the spin splitting of charmonium S-waves does not converge very quickly. [TT 



In the present work, the spin splittings of S-wave heavy-light charmed mesons nicely satisfy 
|0(1/M 3 )| < \0(1/M 2 )\ < \0(1/M)\. However, there is an individual 0(1/M 3 ) term which 
is larger in magnitude than the total 0(1/M 2 ) contribution. Future studies may be able 
to improve this situation, and various suggestions present themselves as conclusions to this 
exploratory study. 



II. ACTION 



The lattice action can be written as the sum of three terms, 
S = S G (U) + S q (q, q; U) + Sq(Q, Q; U) 



(1) 



where U, q and Q are the gauge field, light quark field and heavy quark field, respectively. 

Following the work of Liischer and Weisz fll2| |, a gauge field action which is classically- 
correct up to 0(a 4 ) errors can be written by including a sum over 1x2 rectangular plaquettes 
(Urt) as well as 1 x 1 elementary plaquettes (U p i), 
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A tadpole factor, defined by 
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U = (^ReTrf^) 1 / 4 



(3) 



has been introduced to absorb the lattice tadpole effects and thereby improve the matching 
to perturbation theory ||T3| . 

A light quark action, with classical errors at 0(a 2 ) in spectral quantities, has been 
constructed by Sheikholeslami and Wohlert [ 13 ], 

S q {q,q;U) = -^2q(x)q(x) 
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Again, the tadpole factor has been included for the reduction of quantum discretization 
errors. The lattice field strength tensor is given by 



gF^x) = I (fi^(x) - Q\ v {x)) - ^Im(Tr^(x)) 
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(5) 



U^{x)U u {x + fi)Ul{x + v)Ul(x) 

+ U u (x)Ul(x - fi + £)E^(:e - fi)U^(x - p) 

+ — p)Ul(x — fx — i>)U^(x — fi — v)U v {x — v) 

+Ul(x-v)Up(x-v)U„(x + fi-i>)Ut(xj . 



(6) 



For quarkonium, the form of the heavy quark action has been discussed in detail by 
Lepage et al. @[. It is convenient to write the heavy quark action in terms of the Hamiltonian, 
H, 



S Q (Q,Q;U) = J d A xQ\x)(iD t -H)Q(x) . 



(7) 



To discuss heavy-light mesons, it is appropriate to reorganize the velocity expansion of H, 
discussed in Ref. ||, into an expansion in powers of the inverse heavy quark bare mass, M, 
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The coefficients of the Hamiltonian are chosen so the dimensionless parameters, q, are 
unity at the classical level. Terms arising from quantum effects, i.e. containing powers of 
g unaccompanied by E or B, have not been shown. A tilde on any quantity indicates that 
the leading discretization errors have been removed. In particular, 
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The various lattice derivatives are defined as follows. 
aAiG(x) 
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All of the terms in Eqs. (|9|-[l2D appear in Ref. @ except for the pieces of the eg and C\q 
terms quadratic in B (because they are of negligibly high order for quarkonium) and the C\\ 
term (to be discussed below). The fact that the Hamiltonian H is complete to 0(1/M 3 ) in 
the classical continuum limit has been shown by Manohar [ 15 ]. 

It is conventional to separate H into two pieces, H and 5H, such that the evolution of 
a heavy quark Green's function takes the form P,|T5T| 
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where n is a parameter which should be chosen to stabilize the numerics. Notice that for 
a free heavy quark field, the only relevant terms in the Hamiltonian are H and the terms 
containing c 1; c 5 , c 6 or c\\. Setting each of these to its classical value of unity, and working 
consistently to 0(1/M 3 ), gives 



G 



t+i 



exp 



for t > 



(26) 



which displays the absence of discretization errors for the free quark Hamiltonian up to 
0(a 2 ). The terms containing explicit powers of "a" (05, c$ and en) were added to the 
Hamiltonian precisely for this purpose. 



III. CORRELATION FUNCTIONS 

A heavy-light meson is created by the following operator, 

^Q^x)T(x)q(x) , (27) 

X 

where T(x) is a 4x2 matrix containing the spin structure 

l S : T(x) — (01), (28) 
3 S 1 : T(x) = ( a t ) . (29) 

Gauge-invariant smearing was also tried according to the method described in Ref. ||11|| , 



but it provided no significant improvement for the heavy-light S- waves, which already display 
clear plateaux for local sources and sinks. 

The general form of the meson correlation function is 

G mcson (p,t) = J2 Tt [^(M^y (y - x)^Y\ sk) (y)G t {y - x)T (sc) (x)\ exp (-ip ■ (y- x)) . (30) 

y 

Because NRQCD is an expansion in the inverse bare heavy quark mass, all meson mass 
differences can be obtained from correlation functions at p — 0, but the absolute meson 
mass itself remains undetermined. One way to fix the mass is to compute the change in 
energy when a meson is boosted, 

E >- E ^^i- (31) 

This defines the kinetic mass, M^ n , which is interpreted as the meson's physical mass. For 
the present work, E p is computed only for the 1 S'o state, with p = (0, 0, 2tt/L s ) where L s is 
the spatial extent of the lattice. 
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TABLE I. Simulation parameters. Njj is the number of gauge field configurations, k c is the 
hopping parameter at the critical point, k s is the hopping parameter corresponding to the strange 
quark mass (from mx* /m^=1.8) and a p is the lattice spacing derived from the p meson mass. 



Lattice Njj 


P 


K 


Uo 


K c 




a p [fm) 


8 3 x 14 400 


6.8 


0.135, 0.138, 0.141 


0.854 


0.1458(1) 


0.1398(4) 


0.260(6) 


10 3 x 16 300 


7.0 


0.134, 0.137, 0.140 


0.865 


0.1434(1) 


0.1385(3) 


0.225(8) 



IV. RESULTS 

Gauge field configurations, periodic at all lattice boundaries, were generated using a 
pseudo-heatbath algorithm. After 4000 thermalizing sweeps, the retained configurations 
were separated from one another by 250 sweeps. Light quark matrix inversion was performed 
by a stabilized biconjugate gradient algorithm, also periodic at the lattice boundaries. 

Although the light quark field is periodic in each space-time direction, Eqs. (|24T - |25|) 



indicate that the heavy quark field is periodic only in the spatial directions. Therefore, the 
correlation functions for heavy-light mesons are only useful for times smaller than about 
Lt/2, where L t is the temporal extent of the lattice. Fig. [I| shows examples of effective mass 
plots for 1 S'o and 3 Si charmed mesons, where 

m " (T) = { C—(flr) J ' (32) 

In all plateau is found where the effective mass values computed at three (or more) 

neighboring timeslices are equal within the bootstrap errors. 

The mass of a 1 So heavy-light meson is taken to be the average of all effective mass values 
within the plateau region. The uncertainty associated with this mass is determined by a 
bootstrap procedure, where 1000 bootstrap ensembles are chosen from the original data such 
that each ensemble is the same size as the original data sample. A bootstrap distribution 
is then obtained by computing, for each ensemble, the average of all effective mass values 
within the plateau region. The uncertainty in the meson mass is then taken to be half the 
distance between the 16th and 84th percentiles in the bootstrap distribution. 

The same method could be applied to the mass of the 3 Si, but the errors are expected 
to be highly correlated with those in the 1 S , calculation. Therefore, the 3 Si — 1 S , mass 
difference is determined by applying the bootstrap technique directly to the mass difference. 

Some parameters of the simulations are given in Table |. The determinations of k c , k s 
and a p come from separate simulations involving 250 configurations at /3 = 6.8 and 200 



configurations at /? = 7.0. The values of Uo agree with Ref. [pLlJI . The stabilizing parameter 
n of Eqs. (|2^ - p5|) was at least as large as n = 4 for aM < 1.2, n = 3 for 1.2 < aM < 1.5 
and n = 2 for aM > 1.5. Ol 



The bare charm quark masses at {3 — 6.8 and (3 = 7.0 were obtained in Ref. |jTT| by 
equating the kinetic mass of the r] c with its physical mass. This calculation can be reproduced 
(with poorer statistics) using the gauge field configurations of Table |, provided that one 
additional term is added to the Hamiltonian of Eqs. (|8l-|T2"D, 
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TABLE II. Values for the lattice spacing as derived from the IP-IS mass splitting of charmo- 
nium, and values for the heavy quark masses. Both masses from Ref. [11] lead to M Vc = 2.9(l)GeV. 
All other entries show the range of aM c and aM^ which reproduce the actual experimental values 
for M n , and M r . 



(3 


a hvy [fm] 
(Ref. @) 


aM c aMb 


H sp in (Ref. [11]) H sp in Hf u u 


Hspin Hf u n 


6.8 
7.0 


0.257(9) 
0.205(9) 


1.43 1.5(1) 1.7(1) 
1.10 1.1(1) 1.1(1) 


5.0(2) 5.0(2) 
4.2(1) 4.2(1) 



Hfuii — H — JflQ^fl {^ 2 \ (T • (A x E — E x A) | . (33) 

For heavy-light calculations the term containing c$ is suppressed by four powers of 1/M, but 
in the velocity expansion relevant to quarkonium it contributes at 0(t> 6 ), and is the only 
0(v 6 ) term absent in Eqs. (|]-0). As was done in all other terms, the parameter c$ is here 
set to its classical value of unity. 

It is also worth noting that Hf u u contains terms which are beyond 0(v 6 ) in the velocity 
expansion. The minimal Hamiltonian up to 0(v 6 ) is obtained from Hf u u by omitting the 
portions of the eg and cio terms which involve two powers of B. It is typical to neglect 
the en term as well. In Ref. [ Ill , the entire Cio term was also omitted, since it contains no 
spin structure and therefore seemed negligible for a discussion of quarkonium spin splittings. 
(This point will be addressed below in the context of heavy-light mesons.) The resulting 
Hamiltonian will be referred to as H sp i n - 

Table || shows a good agreement between the bare charm quark masses obtained here and 



those of Ref. [|TTJ. The bare bottom quark mass can be obtained in a similar fashion, using 
the rjb in place of the rj c . Because the rjb has not yet been seen experimentally, it's "physical 
mass" is obtained by subtracting the hyperfine splitting (of the present lattice simulations) 
from the experimental T mass. In practice, the hyperfine splitting is a negligible subtraction 
in comparison with the simulation uncertainties. The resulting values for Mb are shown in 
Table 0. 

Table [H] indicates that for M ~ Mb, the terms which distinguish between H spin and Hf u u 
are negligibly small in comparison with the computational uncertainties. Some evidence of 
their effect might be seen near M = M c , but the large uncertainties do not allow a definitive 
statement to be made. A more precise comparison of H sp i n and Hf u u can be obtained from 
the quarkonium spin splittings, but this is not required for the present study. 

Having fixed all lattice parameters from light-light and heavy-heavy meson observables, 
the heavy-light spectrum will now be considered. Fig. shows the simulation energy of a 
1 S'o heavy-light meson as a function of the bare heavy quark mass, computed to 0(1/M), 
0(1/M 2 ) and 0(1/M 3 ). For each (3, the light quark mass is fixed at a value slightly less 
than twice the strange quark mass, according to Table |. Fig. |^ indicates that terms beyond 
0(1/M) provide small corrections to the leading order result when M ~ Mb, but these 
corrections grow as M decreases. Near M = M c , the effect of 0(1/M 3 ) terms is larger than 
the 0(1/M 2 ) terms. 

In order to understand the origin of such large 0(1/M 3 ) contributions in the charm 
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region, simulations were performed with each 0(1/ M 3 ) term added individually to the lower- 
order Hamiltonian. Results for the simulation energy are shown in Fig. Apparently, all 
terms except the Cio term offer only modest corrections to the lower-order result. In fact, the 
Cio term is unique because it is the only term in the Hamiltonian (up to 0(1/M 3 )) which 
has a nonzero vacuum expectation value. [T7[ This vacuum value is simply an 0(1 /M 3 ) 
shift of the bare heavy quark mass, and therefore produces an 0(1/ 'M 3 ) shift in the meson 
simulation energy. 

The vacuum expectation value can be computed directly from the gauge field configura- 
tions of Table | by averaging over all lattice sites, 

„* / 9 2 rT?2 s 2 A Jl3.6161(20), for (3 = 6.8 

\W 7 " 1 12.0716(16), for (3 = 7.0 . (34) 

If the vacuum expectation value is removed from the action, 



5Sq = St 3 S d4 ^ at ^)§(E 2 + B 2 ) Q6 g 6 (^) 

.2 6 ab I g 2 



— 3-(E 2 + B 2 ) ab - — ( ^(E 2 + B 2 ) 



Q b (x) , (35) 



then the contribution of the cio term to the simulation energy is reduced to the size of the 
other 0(1 /M 3 ) terms, as shown in Fig. f|. Removal of the vacuum expectation value from the 
action is justified because it is simply an addition to the heavy quark mass, which has already 
been removed from the Hamiltonian, Eq. (H), by the standard heavy- field transformation 
of HQET. Mass differences should not depend on whether or not the vacuum expectation 
value remains explicitly in the Hamiltonian. 

Fig. displays the energy splitting between a 1 Sq meson at zero and nonzero 3- momenta. 
The vacuum expectation value has not been removed from the cio term, but its effect should 
cancel in the difference between zero and nonzero 3-momenta. Fig. || shows that the 
0(1/M 2 ) and 0(1/M 3 ) terms offer only small corrections to the leading contribution, so a 
determination of the kinetic mass from Eq. (|31~D will not depend sensitively on the presence 
of these terms. 

Fig. ^ gives the contribution of each 0(1/M 3 ) term to the energy splitting between 
a 1 Sq meson at zero and nonzero 3-momenta. The contribution of each term is small in 
comparison to the statistical uncertainties, including the Ci term. (As expected, both Figs. 
H] and ^ remain essentially unaltered if the vacuum expectation value is subtracted.) 

Fig. |5| raises another important issue. Because the bare quark mass is not a physical 
parameter, it is not necessarily sensible to compare the results of the 0(1/ M), 0(1/ M 2 ) and 
0(1/M 3 ) Hamiltonians at a fixed bare mass. It would be preferable to make the comparison 
at a fixed kinetic mass. However, Fig. [5] demonstrates that the 0(1/ M), 0(1/M 2 ) and 
0(1/M 3 ) Hamiltonians all give rise to the same relationship between the bare and kinetic 
masses, within the statistical uncertainties. 

In Fig. [7|, the spin splitting is plotted as a function of the bare heavy quark mass. The 
vacuum expectation value has not been removed from the Cio term. 0(1/M 2 ) and 0(1/M 3 ) 
terms provide small corrections to the spin splitting near M = Mb, but sizeable ones near 
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M = M c , with an 0(1/M 3 ) contribution that is larger in magnitude than the 0(1/M 2 ) 
contribution. 

According to Fig. |8], there are two 0(1/M 3 ) terms which dominate the large correction 
to the spin splitting: Cj and C\q. The importance of Ciq for the spin splitting is somewhat 
puzzling, since that term in the Hamiltonian is spin-independent. Some insight is gained by 
removing the vacuum expectation value from the cio term ||17|| , which is shown in Fig. ||] to 
remove almost the entire effect of the cio term in the charm region. This small contribution 
of Cio to the spin splitting is to be expected for a spin-independent operator. Apparently 
the large effect of c w in Fig. | is spurious, perhaps because the vacuum expectation value 
is not sufficently small compared to the heavy quark mass. 

To confirm this, recall Eq. (^) which describes the discretization of the heavy quark 
propagator. In the continuum limit, this propagation depends exponentially on the Hamil- 
tonian. At a finite lattice spacing, Eq. (|25|) uses a linear approximation to the exponential 
of aSH. Following Ref. [[|, the exponential of aH was fit more precisely than a simple 
linear approximation, so that an instability could be avoided at small M. In Eq. ([25|), if the 
cio term is subtracted from SH and added to H , then the spin splitting reproduces Fig. |9| 
rather than Fig. [8|. (The values of the parameter n in Eq. fl25|) are not changed from what 
have been used througout this work. For example, n = 3 at (3 = 6.8 and n = 4 at (3 = 7.0 
in the charm region.) This supports the suspicion that the large vacuum expectation value 
was causing a breakdown of the results in Fig. §. The problem has thus been successfully 
overcome in two separate ways: by the explicit removal of the vacuum value from the Hamil- 
tonian, or by a better-than-linear approximation to the exponential dependence of heavy 
quark propagation on the Cio term. 

Fig. | indicates that the spin splitting satisfies \0(l/M 3 )\ < \0{l/M 2 )\ < \0(1/M)\. 
This same ordering persists for all M values considered, as seen in Fig. [Tl| Notice in 
particular that the turnover of the 0(1 /M 3 ) data near M c is completely removed from Fig. 
|7j by correctly accounting for the vacuum value in the Cio term. Meanwhile, the data at 
aM > 2 are not affected in a statistically-significant way. The 1/M expansion might now 
appear to be nicely convergent, but some caution is suggested due to the c 7 term, which is 
by itself larger in magnitude than the total 0(1/M 2 ) contribution to the spin splitting of a 
charmed meson. 

It is interesting to consider the difference between squares of the 3 Si and 1 So masses, 
which empirically is remarkably independent of the "heavy" quark mass: 



(£+*)2 _ (£+)2 ^ ( B 0*y _ (£0)2 = QA8 Ge y2 ( 36 ) 

[B*f - {B s f = 0.51 GeV 2 (37) 

(D + *) 2 - (D + ) 2 w (D *) 2 - (D ) 2 = 0.55 GeV 2 (38) 

(D* s ) 2 - [D s f = 0.59 GeV 2 (39) 

(K + *) 2 - (K + ) 2 w (K *) 2 - (K ) 2 = 0.56 GeV 2 (40) 

( p +)2 _ (7r +)2 w (p o )2 _ (7r o )2 = 5? GeV 2 (41) 



In the extreme heavy quark limit, this result is easily understood using heavy quark sym- 
metry: the spin splitting vanishes as 1/M while the meson masses themselves grow linearly 
with M, so the difference of squares is a constant, 
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m v — m 2 P = {my — mp)(mv + mp) = constant + 0(1/M) . (42) 

For mesons containing only light quarks, the explanation is perhaps not so clear. Some 
authors have related it to chiral symmetry ||18|| . It should be noted that m v —m 2 p rs constant 
is also a consequence of the nonrelativistic quark model with a linear potential (and no heavy 
quark assumptions). 

The difference of squares arising from the present simulations is shown in Fig. |ll]with the 
vacuum value subtracted from the C\q term, and for a particular light quark k. (For all cases 
considered, the results are essentially independent of k, in agreement with Eqs. (|36}{41"|).) 
The large errors are due to the required use of and they are correlated, as evidenced 
by the central values being constant to within a much smaller uncertainty than the quoted 
errors would require. In fact, the 0(1/M), 0(1/M 2 ) and 0(1/M 3 ) data are each constant 
for the full range of M-values that were considered. 



For the purpose of comparison, the same plot is displayed in Fig. but without a 
proper treatment of the vacuum expectation value, i.e. with the vacuum value retained in 
the cio term and the simple heavy quark propagation of Eq. (p5|) . Erroneous results are 
clearly produced at 0(1/M 3 ) for aM < 2. 

For all of the observables under discussion in this work, the en term is essentially irrel- 
evant, which may not be too surprising in light of its large numerical suppression factor in 
the Hamiltonian, Eq. (fT^) . 

Finally, to make the connection to experiment, it is necessary to interpolate to the strange 
quark mass (i.e. to k s ), and to extrapolate to the limit of massless up and down quarks (k c ). 
Interpolations are performed linearly between the two nearest k values, and extrapolations 
are linear in all three available k values. It is also necessary to determine the physical mass 
scale. For (3 = 6.8, both the p meson mass and the charmonium IP — IS mass splitting lead 
to the same physical scale (a p ~ Ohvy)- Use of the bare heavy quark masses from quarkonium 
(Table ||) gives the results of Table 

Atp = 6.8, simulations up to 0(1/M), 0(1/M 2 ) and 0(1/M 3 ) each prod uce masses for 
both the D s and B s which are consistent with the experimental values, indicating that the 
bare charm and bottom masses from quarkonium physics are also relevant to heavy-light 
mesons. The light quark dependence at (3 = 6.8, as probed by D s — D and B s — B, is also 
in reasonable agreement with experiment. The spin splittings are significantly smaller than 
experiment, which is a general feature of previous lattice results as well []§-[|, and is often 
attributed to quenching. 

The situation at (3 — 7.0 is complicated by the fact that a p ^ ah vy . In Table [TT], the 
lattice data are shown for both of these normalizations with the bare heavy quark masses 
fixed to the values obtained from quarkonium. The use of ah vy produces D s and B s masses 
which agree nicely with experiment (as was found for (3 = 6.8), whereas the data normalized 
to a p clearly cannot produce accurate D s and B s masses when the bare masses are fixed by 
quarkonium. 

Conversely, mass differences normalized to a^ vy tend to be larger than the results at 
(3 = 6.8, whereas the a p -normalized mass differences are found to scale remarkably well with 
respect to the (3 = 6.8 results. This preference of the data for a p is in accordance with the 
familiar notion that the dynamics of heavy-light mesons is governed by the light degrees of 
freedom, rather than by explicit heavy quark dynamics. 
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TABLE III. Values for the physical masses and mass differences in MeV. The bare quark 
masses are set to aM c = 1.43, aM b = 5.0 at p = 6.8, and aM c = 1.10, aM b = 4.2 at (3 = 7.0. The 
Hamiltonian contains terms up to 0(1/ M k ), where k = 1,2,3. For charmed mesons, the vacuum 
value has been subtracted from the c\q term. The quoted errors include the uncertainties in k c , k s 
and a. (The uncertainties in k c and n s are negligible except for D s — D and B s — B.) 





D s D a — D D* -D D* s -D s 


B s B s -B B* -B B*-B 




• experiment (Ref. 0) • 
1969 99,104 141,142 144 5369(2) 90(3) 46 47(4) 


k < 1 
k < 2 
k < 3 


• p = 6.8 

2010(130) 96±f 92(6) 84(3) 
2000(110) 99^ 104(7) 95(4) 
2090(120) 98^ 100(6) 89(3) 


• 

5300(700) 75j^ 37(4) 34(2) 
5300(700) 79±^° 39(4) 36(2) 
5300(700) 791?! 39(4) 35(2) 


k < 1 

k < 2 
k < 3 


• (3 = 7.0, scaled 
1940(120) llllii 100(10) 92(5) 
1920(110) 113±X 1 114(8) 105(5) 
2060(110) H3±f 111(7) 100(5) 


by a hvy • 

5000(450) 86l? 40(4) 37(2) 
5000(450) 86^7 43(4) 40(2) 
5000(450) 86±f 42(4) 39(2) 


k < 1 
k < 2 
fc < 3 


• [3 = 7.0, scalec 
1760(100) 1021?! 91(6) 84(4) 
1750(90) 103lg 104(6) 96(6) 
1880(90) 103^ 101(6) 91(4) 


I by o p • 

4500(400) 79±l 36(3) 34(2) 
4500(400) 79±l 39(4) 37(2) 
4500(400) 79±l 39(3) 36(2) 



Perhaps the most satisfactory determination of mass differences at /3 — 7.0 would be 
obtained by normalizing to a p and re-tuning the bare mass to the heavy-light spectrum 
itself, with no reference to quarkonium. However, the possibly-problematic convergence of 
the 1/M expansion for charmed mesons precludes a more detailed effort in this direction at 
present. There is no re-tuning required at (3 = 6.8, so at least in this case an unambiguous 
quantitative comparison to experiment can be made from the data in Table |TJ, although 
concerns about the potentially-large 0(1/M 3 ) contributions (such as the c-j term) must 
certainly be addressed. 

In the work of Ishikawa et al. |J , a range of heavy quark masses were studied and it was 
found that the 0(1/M 2 ) terms increase the spin splitting relative to the 0(1/ M) value, in 
agreement with what is reported here in Table |T|. The classically-complete set of 0(1/M 3 ) 
terms have now been included, and their total contribution is smaller in magnitude than 
0(1/ M 2 ). However, it has also been demonstrated that the small 0(1/ 'M 3 ) contribution 
to the spin splitting results from a cancellation involving an individual term (the c-f term) 
which is dangerously-larger in magnitude than 0(1/M 2 ). 

V. CONCLUSIONS 

The masses of 1 5'o and 3 5*i heavy-light mesons have been obtained from quenched lattice 
NRQCD at two lattice spacings, near 0.22fm and 0.26fm, using tadpole-improved, classically- 
improved light quark and gauge field actions. Results were obtained separately at 0(1/ M), 
0(1/ M 2 ) and 0(1/ M 3 ). The effects of individual terms at 0(1/ M 3 ) were also shown. 
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The simulations up to 0(1/M 2 ) support the existing knowledge of heavy-light S-waves 
for lattice QCD. Masses are in qualitative agreement with experimental data, except that 
the spin splitting is noticeably smaller than experiment. This may be due, at least in part, 
to quenching. 

The contributions of the 0(1/ M 3 ) terms have not been studied in detail previously. A 
novel feature at this order is the existence of a large vacuum expectation value (in the C\o 
term) that shifts the heavy quark mass. Special care is needed when this vacuum value 
is present, particularly in the charm region. An important example is Fig. || where the 
vacuum value seems to make a sizeable contribution to the S-wave spin splitting, but the 
effect was shown to be an artifact of the familiar discretization of exp(Hr) in the heavy 
quark propagation. 

It should be noted that the cj term also produces an 0(1/M 3 ) contribution to the spin 
splitting which is larger in magnitude than the total 0(1/M 2 ) piece. In fact, Fig. |8| indicates 
that this effect is very similar in size to the spurious contribution from the vacuum value 
in the cio term, thus raising the question of whether the Cj contribution might also contain 
some artifact related to the discretization of exp(ifr). This systematic uncertainty has not 
been discussed in the literature to date, but might offer some important insight for the 
charmed spectrum calculation and also for charmonium. 

Other possibilities for reducing the magnitude of the 0(1/ 'M 3 ) corrections also deserve 
further study. In the present work, the coefficients of the NRQCD Hamiltonian have been 
approximated by their classical values, along with tadpole improvement. It would be inter- 
esting to see the effects of retaining one-loop perturbative or nonperturbative renormalization 
for these coefficients. One might also consider working at a smaller lattice spacing, although 
this will move the charm quark mass even further away from the heavy quark limit. (Recall 
that aM c = 1.1 at (3 = 7.0. According to Ref. 0, aM c = 0.81 at (3 = 7.2.) Some benefit 



might come from using a different definition of the tadpole factor, such as the Landau link 
definition, which increases aM c at a fixed lattice spacing. [[H],[BJ 

The existence of an indivdual 0(1/ M 3 ) term whose contribution to the spin splitting 
is larger than the cumulative 0(1/M 2 ) effects, as was found in the present work, indicates 
that the application of lattice NRQCD to charmed mesons requires care. 
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FIG. 1. Effective mass plots for 1 Sq and 3 Si charmed mesons at rest, including terms up to 
0(1/M 3 ). Solid symbols denote data at /3 = 6.8, aM = 1.43 and k = 0.135, while open symbols 
correspond to (3 = 7.0, aM = 1.10 and k = 0.134. The meson source is at r = 1. Horizontal lines 
indicate masses extracted from the plateau regions. 
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FIG. 2. The simulation energy of a ground state heavy-light meson at rest. Results are displayed 
from terms up to 0(l/M k ), with k = 1, 2, 3. M is the bare heavy quark mass. Solid symbols denote 
data at (3 = 6.8 and k = 0.135, while open symbols correspond to (3 = 7.0 and k = 0.134. 
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FIG. 3. The simulation energy of a ground state charmed meson at rest. Results are displayed 
from terms up to 0(l/M k ), with k = 1,2, 3. Solid symbols denote data at f3 = 6.8, k = 0.135 and 
aM = 1.43, while open symbols correspond to (3 = 7.0, k = 0.134 and aM = 1.10. To the right 
of the vertical line, the effect of adding each 0(1/M 3 ) term to the 0(1/ M 2 ) Hamiltonian is shown 
individually. 
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FIG. 4. These data are identical to Fig. ^ except that the vacuum expectation value has here 
been subtracted from the c\q term. 
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FIG. 5. The energy splitting between a ground state heavy-light meson with momentum 
p = (0,0, 2ir/L s ) (where L s is the spatial extent of the lattice) and the same meson at rest. 
Results are displayed from terms up to 0(l/M k ), with k = 1,2,3. M is the bare heavy quark 
mass. Solid symbols denote data at j3 = 6.8 and k = 0.135, while open symbols correspond to 
= 7.0 and k = 0.134. 
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FIG. 6. The energy splitting between a ground state charmed meson with momentum 2tt/L s 
(where L s is the spatial extent of the lattice) and the same meson at rest. Results are displayed 
from terms up to 0(1/M k ), with k = 1, 2, 3. Solid symbols denote data at (3 = 6.8, k = 0.135 and 
aM = 1.43, while open symbols correspond to (3 = 7.0, k = 0.134 and aM = 1.10. To the right 
of the vertical line, the effect of adding each (9(1/M 3 ) term to the 0(l/M 2 ) Hamiltonian is shown 
individually. 
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FIG. 7. The spin splitting of S-wave heavy-light mesons, from terms up to 0(1/M k ), with 
k = 1, 2, 3. M is the bare heavy quark mass. Solid symbols denote data at (3 = 6.8 and k = 0.135, 
while open symbols correspond to /? = 7.0 and « = 0.134. The vacuum expectation value has not 
been removed from the c\q term. 
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FIG. 8. The spin splitting of S-wave charmed mesons. Results are displayed from terms up to 
0(l/M k ), with k = 1,2, 3. Solid symbols denote data at /3 = 6.8, k = 0.135 and aM = 1.43, while 
open symbols correspond to (3 = 7.0, k = 0.134 and aM = 1.10. To the right of the vertical line, 
the effect of adding each 0(1/M 3 ) term to the 0(1/M 2 ) Hamiltonian is shown individually. The 
vacuum expectation value has not been removed from the cio term. 
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FIG. 9. These data are identical to Fig. ^ except that the vacuum expectation value has here 
been subtracted from the c\q term. 
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FIG. 10. These data are identical to Fig. [?] except that the vacuum expectation value has here 
been subtracted from the c\q term. 
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FIG. 11. The difference between squared masses of vector and pseudoscalar heavy-light mesons, 
from terms up to 0(1/M fc ), with k = 1, 2, 3. M is the bare heavy quark mass. Solid symbols denote 
data at (3 = 6.8 and k = 0.135, while open symbols correspond to (5 = 7.0 and k = 0.134. The 
vacuum expectation value has been subtracted from the c\q term. 
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FIG. 12. These data are identical to Fig. 11 except that the vacuum expectation value has 
here not been subtracted from the c\q term. 
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